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We evaluate the three-photon vertex functions at order B and B'^ in a 
weak constant magnetic field at finite temperature and density with on shell 
external lines. Their application to the study of the photon splitting pro- 
cess leads to consider high energy photons whose dispersion relations are not 
changed significantly by the plasma effects. The absorption coefficient is com- 
puted and compared with the perturbative vacuum result. For the values of 
temperature and density of some astrophysical objects with a weak magnetic 

field, the matter effects are negligible. 
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I. INTRODUCTION 



The photon sphtting process in the presence of an external magnetic field was thoroughly 
studied in the vacuum by Adler Since then, it has been studied several times over the last 
years using different techniques and for different regimes [0-0 due to its possible relevance 
in astrophysical applications (see e.g. 0). 

In the vacuum, the Furry theorem forbids the processes with an odd number of photon 
vertices. The lowest order contribution comes from the box graph with an insertion of the 
external magnetic field. In the coUinear aproximation (which is always considered in the rest 
of this paper), this matrix element vanishes due to Lorentz and gauge invariance. Thus, the 
first non-zero contribution to the amplitude comes from the hexagon which is of order B^. 
Adler also studied how the propagation of photons is modified by the presence of a magnetic 
field. He found that there are two linearly polarized modes and due to CP invariance, all 
processes with an odd number of perpendicularly polarized photons vanish. This along with 
the induced birefringence in the polarized vacuum enforces that the only allowed channel is 
||^±_L. The absorption coefficient is 



where 6 is the angle between the magnetic field and the direction of the incident photon, uj 
is the energy of the incident photon and |e|i?cr = ""^^ — 4.4 10^ T. 

In the presence of a hot or dense medium, explicit Lorentz invariance is lost since now 
there is a privileged reference frame, namely the frame in which the system is at rest. 
Therefore the box diagram does not need to vanish. In addition, for a non zero net electron 
density, the ground state is not C invariant and hence the amplitudes with an odd number of 
photon vertices do not need to vanish either. Since the environment in astrophysical objects 
is not exactly the vacuum, it is worth studying the effects that density and temperature can 
have in the photon splitting process. 

As far as we know, a quantitative analysis of the matter effects on the photon splitting 
has been reported in a few papers In Ref. the absorption coefficient is computed 




(1) 
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by including the plasma effects on the dispersion relations but using the vacuum vertices 
from the Euler-Heisenberg effective Lagrangian. In Ref. 0, a one- loop thermal effective 
Lagrangian for a constant electromagnetic field is computed and from that the relevant 
matrix element of order B is read by following the same procedure that is used in the 
absence of a medium In Ref. [g], the same approach is used starting from a two- 

loop effective Lagrangian. Thus, the amplitudes reported in Refs. have an analytic 
dependence on the frequencies of the photons in the process. While this approach works 
properly in the vacuum, it can have a flaw when it is applied for the study of the plasma 
effects. 

In the vacuum, a consequence of Lorentz and gauge invariance is that the effective 
action for the electromagnetic field displays an analytic dependence on the field strength 
and its derivatives. It is just this analytic behaviour which enables a well defined derivative 
expansion of the effective action whose first few terms of the Taylor's series can be computed 
from the Euler-Heisenberg Lagrangian Q Moreover, because of Adler's theorem, the three- 
point photon vertex in the coUinear approximation computed from the Euler-Heisenberg 
Lagrangian is exact to order B^. 

On the other hand, in a medium, Lorentz invariance is lost and the effective action 



displays a non-analytic behaviour |llli|12| . This invalidates the identification of proper photon 
vertices based on an effective Lagrangian computed by a derivative expansion technique. 
In order to clarify this essential point, let us consider the thermal two-point polarization 
tensor in QED at high temperature. In the static limit, the only non-zero component of 
this tensor is H°°(c<j = 0,k) = — -|- 0{k'^), where m^, = is the Debye mass of 

static screening. Then, the corresponding lowest order terms of a derivative expansion in 



^When the masses are different from zero and n > 3, the one-particle irreducible n-point functions 
are analytic functions of all its arguments, pf = tof—pf and pi-pj = uJiUJj—pi-pj for {pi, . . . ,Pn-i) = 
(0, . . . , 0), so that its Taylor series can be guessed from the behaviour for uJi = 0. 
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the effective action are T ~ / d^x [^mlAl + C(E2) + C(B2)) . Here Aq, E and B represent 
the temporal component of the gauge field, the electric and the magnetic field respectively. 
However, when cu ^ k, the non-zero components are ir^(uj = A;, k) = m1{5^^ — A;* A;-'), where 
= is the transverse photon mass defined below in Eq. (^). Now, the corresponding 

effective action is 

r ~ y c^^x ^ [k^ + divA ^ divA) , (2) 

which has a non-local dependence in the fields but is gauge invariant. This differs from 
the corresponding effective action based on the derivative expansion and would not have 
been obtained from an expansion about a constant field. Certainly, in order to describe a 
dynamical process involving propagation of photons with a dispersion relation uo"^ k"^^ the 
relevant matrix elements should be computed from an effective action similar to the one in 
the Eq. (1). 

In this paper, we will directly compute the three point amplitude using the Schwinger 
propagator for an electron in an external magnetic field when the collinear external momenta 
are on shell. Expanding the three-photon vertex at lowest order in the external field and 
retaining the relevant pieces we will obtain the box contribution. Proceeding in the same way, 
we will also obtain the order contribution which would come from the pentagon graph. 
This leads to Eqs. (^) and (^) below, which are the main findings of this paper. Both 
of them display a non-analytic dependence on the photon energies. By computing several 
components of the amplitudes, we have verified that the Ward identities are satisfied. 

The presence of a medium not only affects the amplitudes, allowing some of them not 
to vanish, but also modifies the photon propagation properties. However, at high energy we 
expect this modification to be negligible compared with that due to the magnetic field. In 
order to estimate under which conditions this is true, we recall that the dispersion relation 



of high energy photons in the absence of an external field is a; ~ yP^ + , where rrit is 
known as the transverse mass of the photon |13 
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— / dp—{nF{ep-fi) + nFiep + n)), (3) 
vr JO 



with npix) = 1/ (e^^ + 1) and Ep = y/p'^ + the energy of an electron of momentum p. The 
transverse mass is of the same order as the plasma frequency. The refractive index is found 

2 

to be n 1 — On the other hand, the refractive indices due to vacuum birefringence 
are i 

2 

Brr j ' 907r' " QOtt 



I ^ I / 5 sin ^ . ^ I , , ^ . 

riu = 1 + aa,i — , a|| = , ct = • (4j 



In this paper we follow the convention of Adler |jT|, where -L (||) means the magnetic field 
of the photon is perpendicular (parallel) to the plane containing the external magnetic field 
and the photon's momentum vector. Thus, we must require 

" ~ TSs (bI^) ^ p»io=(%^)'".--o. (5) 

This inequality specifies the condition that is necessary in order to neglect the medium 
effects on the photon propagation. We will find that in these conditions the contribution of 
plasma effects to the absorption coefficient is smaller than the vacuum contribution. 

In the next section we describe the computation of the box and pentagon amplitudes. 
The following section is devoted to obtaining the absorption coefficient and evaluating it in 
several situations. Finally in section IV we give a brief summary. 

II. THE BOX AND PENTAGON AMPLITUDES 

In order to compute the three photon amplitude we need the fermion propagator in the 
presence of an external magnetic field. It was computed by Schwinger a long time ago 
using the proper time method 

G(X, X') = -im -e4 + m) ^ ds , ~l , ^i{m^-ie)s^H-,^') ^i%o>^-F,^s /g^ 

^ ' ^ J~oo (47r)2s2 ^ ^ 

/I i 
X exp I — Tr ln[(eFs)~"'^ sinh(eFs)] H — (x — a;')e-F coth(eFs)(a; — x'] 
\ 2 4 

with 



^x,x')=-ie r dx^{A^{x) + ]-F^,{x-xr). (7) 

Jx' / 

It can be expressed as 

G(a;,x') =e*("'"') ^(x-a;'), (8) 

where G{x — x') is explicitly invariant under translations while the phase factor is not. We 
need the expansion of G in powers of the magnetic field, for which we obtain in Fourier 
space 

Gip) = Goip) + G^ip) + G2ip) + 0(3'), (9) 
ml+ 



Goip) 
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~ P ~ 771 



with F^i, = e^j,Q,^F°^/2 and a^'^ = |[7^,7'^]- We have arrived at these results by expanding 
the integrand of G in powers of B, performing the Fourier transform and finally doing the 
integral over the proper-time parameter. The expansion of G in powers of B had been 



carried out previously in Ref. |T5[, where it was pointed out that the phase factor gives a 



non trivial contribution to the three-point amplitude 

. ^Mx,z) . ^Mz,y) = exp (^-^{z - yTF^p{y - xf ] (10) 
Mx,z) . ^Hz,y) . ^Hy,-) ^ exp f ^(y - a;)"F,^(x - 



exp \^—{x - zYF^p{z - yy 

Any of the three forms can be used in the computation of the amplitude due to the trace over 
the Dirac indices. Since the phase is already linear in B, there are two kinds of contributions 
at first order: one when there are two factors Gq and one Gi and the other when there are 
three factors Go and the phase. 

We start with the contribution without the phase which can be written directly in Fourier 
space (see Fig. 1) 



Mr[p,Pi,P2] =-«Vr^ J ^Jr[YG{k + p,)YG{k-p,)rG{k)] (11) 

To proceed we choose the momentum of the incident photon along the OZ axis: = 
{u!,0,0,p = u) and the magnetic field in the plane OYZ: B = {0,By,Bz) without loss of 
generality. Thus, in the coUinear aproximation the momentum of the outgoing photons are 
Pi = {ui,0,0,pi = ui) and P2 = (072,0,0,^2 = ^^2)- The polarization states are chosen 
following the convention of Adler 0: ef = (0, 1, 0, 0), = (0, 0, 1, 0). 

Introducing the expression for G{p) in Eq. ( [TTl) and taking the linear part in B we can 
compute the amplitude for the different channels. We have written a Mathematica code 
to carry out the Matsubara sums and the angular integrals. The Matsubara sums lead to 
derivatives of the Fermi distribution function which are integrated by parts, resulting at the 
end in the single integral of Eq. (|^) defining the transverse mass. We obtain 

Ar""=-i^-™f. (12) 



M-"^ = i4^m?. (13) 



Mt ^" = -^-m2, 14 
2 a;2 



1 eBy e^ujf + ijj\ijj2 + ijjl^ 
Mj'^" ^ = Mj'^^" = Mi^^" " = M^^^^ = 0. (16) 



/ UJUJ1UJ2 



Now, in order to compute the contribution of the phase, we write the three-point function 
with the external legs in position space 



r^^^[x, y, z] = -{T^'e' J d^x'dyd^z'Ti[Y'Goiy' - x')Y'Go{x' - z')Y'Goiz' - y')] (17) 
xD,y (x - x')D^,^iy - y')D,y{z - z') i^-^{z' - yTFaf^iv' - a^')^) + ( /i, 1/ ^ z^, ^), 

where Da,a' {x — x') are the propagators of the photons. In Fourier space this is given by 

T'^^-^ix, y, z]= J (21^4 ^2!!^ (l7^ ^^.P'(P)^M,.'(Pi)^...'(P2)e-^-e--^-^e--^^- 

xi27r)^5\p - p, - P2)MP'^'''' {p,p^,P2). (18) 



Thus, the corresponding three-point amphtude reads 



(f'k ie d d 



Mr\v.VuV.\ =-(0^e3rE / ^'^f'' (19) 



Proceeding as we did before for the other contribution we obtain an identical result. Thus 
the total amplitudes are 

M|Hxx = -^-m?, (20) 

MxHI^ = ^-^?> (21) 
uji 

M^_x|| = ^-m?, (22) 

eBy e{uj\ + uj^U2 + ^j) ^2 /9o^ 

M||^ll s_ = M||^x|i = Mx^ii II = Mx^x± = 0. (24) 

The amplitudes with an odd number of perpendicularly polarized photons vanish. These 
results can be written in a compact form 



2mf [ 1 1 1 



^IJ yk ^jk yi ^ ^ki yj 

UO2 OJ UOi 



(25) 



where = e^^^BkPm and i, j, k can take the values 1, 2 for polarizations ||, ± respectively. 

Now we turn to the contribution at order coming from the pentagon graph. The 
computation proceeds along the lines of the previous one. The contribution at second order 
in B with no phases is computed as before. For the term linear in the phase we start from 
Eg . (p!9D , but now the factor inside the trace must be replaced by its first order expression in 
B. Finally, for the second order contribution of the phase we have an expression similar to 
Eq. (pD 

Tr[(7^Go(A; + pi)7^Go(A;-j92)7'^Go(A;))] + ^ u.p^). 

The results of each part alone are not very illustrative, so we give only the total amplitude 
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M,HU = -.e^(^)(-;^jn.. (27) 
Af,^,,=-,e»f^Ui*;^V„„ (28) 



with 



. 3 (ByBA (ujf+LJiLJ2 + UjjY . . 

M±^±± = -le rie, (30) 

M||^xx = Mx-.||± = M^^xii = M||^j| II = 0, (31) 



ne = 2 I J2^^'^^^^P - f^) - ^pi^p + /u)) (32) 



being the net electron density. The amphtudes for an even number of perpendicularly 
polarized photons vanish. These results can also be written in a compact form 

M^^\p,p,,P2) = -le' ne (33) 



cr 



X 



U;^ UJ1UJ2J VCoif 00002 1 V^2 OOUOiJ 



where B\ = B^ — (B ■ p)p* and i,j, k can take the values 1, 2 for ||, ± respectively. 

An important issue concerning the outlined computations is its consistency with gauge 
invariance. Following the same procedure, we have computed other components of the three- 
point amplitude and we have explicitely verified that the Ward identities in the collinear 
approximation for null vectors, M^^*^ — M^'^'^ = 0, are satisfied. 



III. THE ABSORPTION COEFFICIENT 

In order to obtain the absorption coefficient a phase space integration with the proper 
measure has to be performed. When the effects of the magnetic field are considered in 
the dispersion relations, Adler showed by using arguments based on CP invariance and 
kinematics, that in the vacuum the only allowed channel is ||— >-L -L. As we have mentioned 
above, for high energy photons {oo ^ loq) propagating in a magnetized plasma, the infiuence 

9 



of density or temperature on the dispersion relation is small and we can expect that the same 
channel is the only one allowed for contributions which are even in the chemical potential. 
For this decay we have 

where the enhancement Bose factors appearing in the probability have been replaced by 
one, because we consider the case uji,uj2 ^ T. If dispersion due to the magnetic field is 
also neglected, the integrals reduce to dki, which can be approximated by the overall 

integral, so that 

a f B sin 6 \^ f f mt\^ 



- ^- (35) 



327r LJ 8 \ Bcr J \Ld J \m J 
For the pentagon contribution the kinematical constraints are the same as before, but the 
CP arguments now allow only procesess with an odd number of perpendicularly polarized 
photons. When both constraints are taken into account the only allowed channel is ||— s>||i±2 
or II II 2. The absorption coefficient is given by an expression similar to Eq. (0) with 
|^IH±±P replaced by |M||^||^p + |M||^_|_||p. Now, the integrand depends on ui and UJ2 and 
hence the final result displays an explicit dependence on the cutoff 



sin=(2») (^] (r^X m ra. (36) 



3 \BcrJ \ujj \uJoJ , 

We stress again that an application of the above formulae for low energy photons is mislead- 
ing, as the formulae themselves show. At low energy, the propagating modes in the plasma 
are longitudinal or transverse plasmons whose dispersion relations differ significantly from 
that we have assumed, u = p. 

We show the explicit dependence on the plasma parameters for two opposite regimes. 
For a completely degenerated electron plasma we have 



a I B sin ^ \ / m 



In 



m I m\ m? 



m, (37) 



3 / R \ / /m \ 2 / mn \ 3 / / « 2 rrr-t 2 \ 3 
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and for high temperature T ^ m and /i = 0, we obtain at order 

Obviously, there is no contribution at order B^. 

In order to find an estimate of matter effects in photon sphtting, we compute the absorp- 
tion coefficient for two astrophysical objects. For a white dwarf we consider the following 
values of the parameters flBl 



— ~2■10-^ — ~2■10-^ ^- 1.02. (40) 

Bcr m m 



The transverse mass and the charge density can be estimated with formulas ^2]), giving 
~ 5 ■ lO^^m, He ~ 2.7 • lO^'^m^. In order to satisfy (|]) we need u ^ 2.5 ■ 10^ m. Taking 
6 = 7t/A and u ~ 10^ m we get the following estimates 

~ 10-^2 m, ^ 10"^^m. (41) 

The value in the vacuum from the hexagon is k\i^±± ~ 3.5 ■ 10^^^ m. Here, the effects of 
matter are completely irrelevant. For a neutron star we can take the values p!6[ 



B T LL 

2-10-\ 1, -^~ 600. (42) 

Brr rn m 



Now we have rrit ~ 40.9 m, ~ 7.3 ■ 10^ and c^q ~ 2 ■ 10^ m. Again, taking 6 = ir/A and 
uj ~ 10^ m we get the estimate 



- 5- 10"'^m, ~ 1.4-10"'^m, (43) 

while the value of the hexagon in the vacuum is k\\^±± — 3.5 ■ 10^ m. We see that the effects 
of the medium are neglible by many orders of magnitude. In addition, the results for the 
neutron star do not change appreciably in the range of temperatures from T/m = 10~^ to 
T/m = 1. 
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IV. SUMMARY 



In this work we have computed the contribution of order and 5^ to the photon 
sphtting process at finite temperature or density in the colhnear aproximation. We have 
computed directly the three-photon amphtude using fermionic propagators in the presence 
of a magnetic external field for photons with uj = p. The phase of the propagator contributes 
non trivially to the amplitude. In order to have photons with ~ p we have restricted our- 
selves to photons with uj ^ ujq so the effects of the medium on the dispersion relations can 
be neglected compared to those of the magnetic field. This condition leads to consider very 
high energy photons, for which the effects of the medium are neglible compared with the 
perturbative vacuum contribution coming from the hexagon, as seen in explicit computa- 
tions. It is not excluded that for lower energy photons, when the full effects of the dispersion 
relations must be taken into account, the box contribution is more significant. 

A natural extension of this work is the study of the matter-induced three-point vertex at 
strong field, a regime with possible astrophysical applications. At strong field, the modifi- 
cation of the dispersion relations due to vacuum effects play a more significant role than 
at weak fields. 
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FIG. 1. The photon sphtting amphtude. The sohd lines represents fermionic propagators in the 




presence of a magnetic external field. 
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